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We model the effects of a large number of zero modes for Nf species of quarks at finite vacuum 
angle 6, using a matrix model with gaussian weights constrained by the topological susceptibility 
and compressibility. The quenched free energy exhibits a cusp at 6 < n that is sensitive to the 
accuracy of the numerical analysis and the maximum density of winding modes. Our results bear 
0^ , much in common with recent lattice simulations by Schierholtz and others. The unquenched free 

energy exhibits similar sensitivities, but for small quark masses or a large density of zero modes 
the results are in agreement with those derived using anomalous Ward identitities and effective 
Lagrangians. 

o : 

i . 

I. INTRODUCTION 

m ; 

QCD with a finite vacuum angle is subtle. Canonical quantization and variational calculations || suggest 
that the vacuum state depends on 9, while covariant quantization seems to indicate otherwise The issue of the 9 
, angle in QCD and the U(l) problem are intertwined [§j4[|^]. At finite 9, QCD breaks CP. Due to the U(l) anomaly, 
the 9 term may be traded from the gauge fields to the quark mass matrix. Bounds from the neutron electric-dipole 
O I moment yield 9 < 1CT 9 [§. 

0^ ' In QCD the dependence of the vacuum partition function on the 9 angle involves an understanding of the vacuum 
0^ ' physics which is essentially nonperturbative. First principle calculations are limited and difficult. The reason is that 
(-h , most lattice QCD simulations rely on important samplings by Monte Carlo techniques which require positivity of the 
^ action configuration by configuration. At finite 9 the action is complex in Euclidean space. 

Lattice simulations using CP n models as well as Yang-Mills theories [Q-^] have been recently carried out, with 
somehow oppos ite conclusions. A number of effective models have been used to gain insights to this important 
problem ||ni|~|l2||, including recent conjectures ||l3|| . Unfortunately the conventional lore of power counting, such as 
chiral perturbation theory, becomes subtle at finite 9 J^Jl^l , although some exact constraints can still be inferred from 
Ward identities @,?,|. 

In this paper, we will assume that the vacuum supports a 9 angle and proceed to analyze some related issues using 
a matrix model. The main thrust of our investigation is to try to find out the conditions under which a numerical 
analysis of this problem (albeit in a model) compares to analytical or quasi-analytical solutions. In sections 2 and 3 we 
introduce the model, and discuss the quenched case. The free energy is found to depend sensitively on the maximum 
density of winding modes [] n and the accuracy of the numerical calculations. Our observations are similar (although 
not identical) to those reached recently by Schierholtz and others using lattice simulations. In section 4, we 

discuss the unquenched case and show that under general conditions the saddle-point results agree with the numerical 
calculations. Our conclusions are in section 5. Some technical details are given in the Appendices. 

II. MATRIX MODEL 

Consider the partition function described by 



1 In the unquenched case the number of winding modes is the number of zero modes. 
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^)-(n«(^""^)>- 

where the averaging is carried using the weight 

f dWdWh-i nTTW " w e^e^ . (2) 

Here W is a complex hermitian asymmetric n + x n_ matrix, n — n + + rt_, and a±x — 2n± — (n). The mean number 
of zero modes (n) is either fixed from the outside or evaluated using the gaussian measure (|J). For simplicity in this 
paper we used the quenched measure without the fermion determinant to fix (n). Throughout, the value of the quark 
condensate E = 1 in the chiral limit. It is readily reinstated by dimensional inspection. 

For a recent review on matrix models in QCD we refer to Ga] (and references therein). In short, (0) with gaussian 
weights is borrowed from the effective instanton vacuum analysis JToj ] where n + counts the number of right-handed 
zero modes , and n_ the number of left-handed zero modes (restricted to zero dimension). The number of exact 
topological zero modes is commensurate with the net winding number carried by the instantons and antiinstantons. 
By analogy with ||l6| , x* an d <J* will refer to the unquenched topological susceptibility and particle compressibility, 
respectively. Higher susceptibilities may be enforced by non-gaussian weig hts provided that the QCD beta 

function is restricted to its one-loop form. This alternative will be discussed elsewhere. 

From [|l6| we have a 2 = 4n„/6 with b = UN c /3 and n» = (n)/V the mean density of zero modes. If the 
compressibility er* is assumed small in units of E = 1, then typically n ~ (n). We note that for the canonical choice 
n„ = 1, <r„ = 0.6 for N c = 3. Finally, through a chiral rotation the 9 angle may be removed from the determinant to 
generate an extra phase e l%s in the measure (see also (pT|)). In this form, the 27r periodicity in 9 is manifest. We recall 
that in the original instanton model, CP is explicitly upset at finite 9, with the exception of 9 = ±n (mod 2tt) 

III. N F = 

In this section we will only analyze the quenched partition function with Nf = 0, thereby probing the nature of 
the measure (^|). As we will show, this is not a trivial exercise and the outcome bears much in common with current 
quenched lattice simulations. First, we discuss the case where (n) = oo with no restriction on the value of n, and 
hence no restriction on the value of x- Second, we discuss the case where (n) is large but finite, so that |x| < n with 
typically n ~ (n) for a peaked distribution in n. 

A. Infinite Sum 

When the sum is unrestricted and infinite, we have for the quenched partition function (up to an irrelevant nor- 
malization) 

oo 

Z Q (6) = Z(0,O) = e lx9 e-* 2/2Vx * . (3) 

X=-oo 

Using Poisson resummation formula we have 

Z Q (0)= J2 e-i v ^ 9 - 2 ^ 2 =6 s (6/2,e- T ). (4) 

k— — oc 

where the last equality involves the third elliptic ^-function with r = 1/(2V%*). The result is manifestly 2n periodic. 
The vacuum free energy, Fq{9) = ~IhZq(9)/V as V — > oo is simply 



2 Under CP 9 = -it goes to — -it which is the same as ir because of the 2tv periodicity. Mutiple degeneracy of states can however 
take place 
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F Q {6) = min-x* (0 + mod 2tt) 2 (5) 

in agreement with the saddle-point approximation to ( H) This simple result is the same as the one obtained using 
large N c arguments prf| , and recent duality arguments |13|. We observe that the cusp at 9 = ir (mod 2tt) sets in for 
V = oo. For finite V the sums converge uniformly, so that 

F Q {6) = -yi^f)^ Xsin( X 6)e-^/^. (6) 

which is always zero at 9 = tt. As V — > oo, the interchange of the derivative with the sum is not valid, hence the 
cusps. When translated to the quenched instanton calculations, these cusps indicate a new phase with spontaneous 
CP violation. A similar observation was made in the context of 1+1 compact electromagnetism using the character 
expansion p9| . 

Finally, we note that (0) can also be rewritten as 



-, I — / +oc +oo \ 

Z Q {9) = Z E {8) + Z {9) = \J*-[ £ e-^- k ^+ £ (-D* e -£(*-M 2 

\k— — oo k— — oo / 



e-^)+9 2 {9,e-' lT ) (7) 



where Ze is n periodic and Zo is n antiperiodic. The appearance of (— l) fc in the second sum in (Q) is important for 
restoring the 2tt periodicity in the full sum. In the thermodynamical limit, the even and odd sums are dominated by 
single Gaussians with F E = x *{6- kn) 2 /2 for \9/tt - k\ < 1/2, and F = F E + ikn/V. 

B. Finite Sum 

In the model we are considering the sum over x is restricted to |x| < N, with A = max n. We will denote by 
n = N/V the maximum density of winding modes. While in general n ^ n*, for a peaked distribution in n (small 
compressibility cr* ) we expect n ~ . This will be assumed throughout unless indicated otherwise. Hence 

N-i 

Z Q {9)= e iex e-* 2 / 2Vx * (8) 

X =-(iV-l) 

Approximating the sum in (||) by an integral and evaluating it by saddle point we obtain Zq ~ e ~ v x*s /2 _ 

In Fig. we show the numerically generated result versus the saddle point approximation for the full free energy 
(27r-periodic) and different values of n. The normalization was chosen so that -Fq(O) = 0. For N = 250 and x* = 1> 
the double precision (16 digit) numerics (circles) breaks away from the saddle point approximation (solid line) at 
9/n ~ 0.2 for both n = 1 and n = 4, while the high-precision |2(J (64 digits) calculations (dashed line) break away at 
9/n ~ 0.3 at n = 1 but agree with the saddle point result at n = 4. 

The leveling for small values of n persists even at infinite accuracy, and in our case is caused by the finite range of 
the summation over x m (§)• Indeed, using the Euler-MacLaurin summation formula, we have 

n-i r N ^ i i 

E f^) J o f(x)d X = - 5 [/(0) + f(N)} + -[/'(AT) - f(Q)] ^[f"(N) /'"(0)] + . . . (9) 

where . . . stand for odd derivatives of /(A) = 2 cos(N9) e - N2 / 2V ** . The generic behavior of the correction terms is a 
non-exponential prefactor times e~~ N / 2 ** v with no dependence on 8 |J This should be compared to the leading order 
result, e~ Vx * e / 2 , hence, a breakdown of the saddle-point approximation for A = V is expected at x*® 2 ~ Vx*j that 
is 9/n ~ l/x*7r « 0.3, in agreement with the high-precision calculations of Fig. |l|. 

If we were to increase the density of winding modes to n > nx*, then for say n = 4, the saddle-point result is 
recovered as indicated by the stars (Fig. [l]). The breaking point noticed above, now lies outside the period of the 
free energy (at 9 ~ n/x*)- Again, the precision in the numerical calculation is important. For n — 4 the precision 
is upgraded from 100 digits for 9/n < 0.3 to 150 digits for 9/tt > 0.3. We have observed that for low precision 
measurements (16 digits), the numerical calculations deviate from the saddle-point approximation for small values of 
9/tt r~j 0.2 (circles) even for n = 4. 



3 8 enters only in the prefactor which drops when taking the logarithm and dividing out by V. 
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C. Comparison to Lattice Simulations 



The observations of the preceding paragraph may be summarized as follows: the leveling of the quenched free energy 
of the matrix model as a function of 9 is sensitive to the numerical accuracy of the calculation. The leveling stabilizes 
at large numerical accuracy, and is found to depend on the maximum density of winding modes n (typically n ~ n») 
and its relative magnitude to the topological susceptibility. For the measure (^), the leveling occurs at 9* ~ n/x* in 
agreement with analytical estimates. 

In an interesting series of investigations, Schierholtz Hj analyzed numerically the effects of a finite 9 angle using a 
CP 3 model in two-dimensions, and also Yang-Mills theory in four-dimensions. The CP 3 simulations shown in Fig. || 
(left) indicate a leveling of the free-energy for 9 < tt/2, cautiously interpreted as a possible evidence for a first order 
transition to a CP-symmetric state Q. We note the striking similarity of this Figure with Fig. [I]. 

These findings were recently reexamined by Plefka and Samuel S who concluded that the apparent first order 
transition was a possible artifact of the accuracy of the numerical simulation for fixed lattice size, as shown in Fig. |^ 
(right). Our findings in the matrix model confirm this observation. The larger the volume V, the more precision 
is needed (exponential precision for infinite V), due to large cancellations in the partition function caused by the 
oscillating phase e lxS . However, we have also found that the leveling depends quantitativaly on the maximum density 
of winding modes n and persists whatever the precision for n ~ 1. In this sense, it would be very useful to understand 
the dependence of the results in on n, with in particular the ones shown in Fig. @. 



IV. N F > 



To assess the effects of light quarks on the partition function of the matrix model, we will consider in this section 
the general case with Nf > 0. We will evaluate the free energy of the matrix model in the saddle-point approximation 
after bozonization, and compare the outcome to direct numerical calculations using large ensemble of asymmetric 
matrices or quasi-analytical methods. 



A. Bosonization 



In (0) the fermion determinant can be rewritten as an integral over A/n-component Grassmanians tp = 
with dimensions {n + ,nJ), that is 



dot j 



(10) 



The 9 angle can then be removed from the action by a U(l) transformation tp — > e nse ^ 2NF q with 75 = diag(l„ + , — l n _). 
The Jacobian of this transformation is just 

det F e^ e ' N f = e «>(»+-™-) (11) 
which is the well known trade-off through the U(l) anomaly. As a result, we have the following Ward identity 

(n+ - n_) B = m (q i H 5 q) g ■ (12) 

Assuming that the range of resummation over \ = n + ~~ n - is infinite, that the distribution in n = n + + ri_ is peaked 
(n — » (n)) and trading the x _ sum by an integral in (uh, we have, 



Z(9,N f ) = / dPdP^expl (n) 



-Tr|P| 2 + iTrlog|z + P| 2 



2n* 



log det— 



z + P 



+ Pt 



(13) 



where z = diag nije t0 / N f . The logarithm in ( |l3| ) is multivalued and its 2ir determination will be assumed. In analogy 
to the quenched case, the approximation of the sum by an integral should be valid for n ~ n* > trxtop, where we 
substituted the quenched topological susceptibiliy x*i by the unquenched one Xtop < X*- F° r small masses Xtop gets 
screened (see Appendix C). 

To keep a tab on the validity of trading the x-sum by an integral, we observe that for Nf = 1 and z — me' 9 , an 
exact form can be reached for 
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Z{0,l) = ({z + P) n +{z + P^) n -) (14) 

followed by the substitution P — ► W in the measure @j. Here P is a complex variable. Because of (|ll]), the same 
partition function can be written as 

oo 

Z(6,l)= e lxe e~ x2/2Vx *Z x (15) 

X=-°o 

where Z x is the partition function for fixed asymmetry x- Again for a peaked distribution N ~ (n), so that 

N- x + 2 p / _/V + x + 2\ 

^=(|) 2 ^e-^A_^ lFl( ^±|^ ;X+1;i v m V2). (16) 

where \F\ is Kummer's (confluent hypergeometric) function. For even x, it reduces to an associated Laguerre 
polynomial in Nm 2 . We note that for N — > oo with Nm and x fixed, Z x ~ I x (Nm) which is the expected generating 
function for the microscopic sum rules plj . 



B. Saddle-Point Approximation 



Without loss of generality, we can set P = diag pje^ e ^ Nf , so that the unsubtracted free energy associated to 
@ reads 



n Nf 

3=1 



X* 
2 



arg 



n 



v 



Pje 



-i<)>j 



\ 



3d2?r 



We note that for Nf — ( |l7| ) reduces to (JsJ) . For small masses rrij <C x* ~ Pj 1 ^ 1 an d (0) s i m plifi es to 



(17) 



3=1 



cos 0j + -x* 



^ + mod 2tt + 51 — 



sin ( 



3=1 



\ PJ 



(18) 



to order 0(m 2 ). The saddle point in the p's decouples and gives pj = 1 — m,/2cos0 + 0(m 2 ). To the same order, 
the saddle point in the 0's is 



= 5^ (^j + 0(771) and TOisin^i = ... = mishit 



>JV> 



(19) 



3=1 



These equations were derived using large N c arguments (T^JT^] and anomalous Ward-identities 
For Nf = 1, we have 



X* 



for which the subtracted free energy is AF = n„m(l — cos( 
For Nf = 2, we have 



1 sin6» + 0(TO 2 ) 



(20) 



sin^i 2 = ± 



7712 isinf 



\J m\ + m 2 + 277117712COS 9 



and the subtracted free energy now reads 

1 



11* 



AF(9) = \nii + m 2 \ — yj m\ + m 2 + 2mi77i2 cos 6* . 



(21) 



(22) 
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For mi = m,2 a cusp develops at 9 = it, since AF(9) = n*|m|(l — |cos#/2|). In this case both numerator and 
denominator vanish in (^l[), hence any value of <f>\ and 02 is allowed provided that 0i + 02 = T- A similar behavior 
was noted by many fl5T ]ll|]l2| , |2"2] , p3] ] , follo wing the spontaneous breaking of strong CP. 

For Nf — 3, the explicit solutions to (|19| ) are in general involved, thereby making an analytical form for the free 
energy involved. At 9 = n, however, the analysis simplifies. Using ( |l9| ) we obtain a trivial solution with one of the 
0's being tt and the others zero, and a non-trivial one, 

cos 03 = 3 13 2 ^ (23) 

2TO3TOim2 

which is doubly degenerate for 

7711771.2 > TO3ITO1 — 7772I . (24) 

The trivial solution corresponds to no cusp at 9 = ir, while the non-trivial one yields a cusp at 9 — -k because of the 
double degeneracy. Again similar observations were made using effective Lagrangians |Il"| , p2| , p3l a nd anomalous Ward 
identities H , where J24|) is known as Dashen's condition p£| . For sufficient flavor breaking (p4j) is not fulfilled and 
strong CP is not spontaneously broken. This is the case in nature where 7773 3> mi, 7772. 

C. Numerical Analysis for Nf = 1 

In Fig. H (left) we compare the fixed size results for the free energy with (77) = 14 (circle) and (77) = 100 (plus) 
over one period Q using single precision (16 digits) numerical calculations for 777 = 0.5, x* = 1. The saddle point 
solution is indicated by the solid line. Increasing the size of the matrices worsen the agreement with the saddle-point 
result (solid line). The reason is that for 7r/2 the integrand in the sum develops alternating signs with large numerical 
cancellations. The breakdown of the numerical calculation takes place at ~ — (In e)/((n)x*) with e the numerical 
precision. These observations are new as they pertain to the unquenched free energy. 

In Fig. U (right) we show the same results after averaging over the size of the matrices (with er*=l), restoring the 
2ir periodicity. Changing the gaussian measure in 71 to a uniform one with \n — (77) | < 6 does not make a noticeable 
change on the Figure. We compare the low (16 digits, circles) and high (64 digits, boxes) precision calculations for an 
averaged size (77) = 14 using an ensemble of 500000 matrices. Large size samplings at large 9 requires exponentially 
larger precision. Much like the quenched case we observe the same dependence on n ~ n„ causing a departure from 
the saddle point results at 9 ~ 0.67T for n = 1. 

Since the n > 1 region is unaccessible by direct use of random matrices, we will use the partition function for fixed 
asymmetry Z x in ( fL5| ) in the saddle point approximation and perform the final sum in ( |l5| ) numerically. We have 
checked that the saddle-point result agrees with the result obtained by averaging over large ensembles of matrices 
with fixed x for n — 1- We refer to this analysis as quasi-analytical and the results are shown in Fig. || (right) for 
n = 1 (dotted line) and n = 4 (dashed line). 

For small masses and 9 < tt/2 the free energy in the saddle-point approximation is F(9) ~ n„m(l — cos#) with no 
cusp whatever 9. For large m & cusp at 9 — n develops, following the decoupling of the flavor (roughly 777 ~ NfX*)- 
In Fig. I the occurence of a cusp for m = 0.5 and n = 4 is an artifact of the quasi-analytical analysis where only 
the saddle-point in x is retained (see above). Indeed, the numerical results from this procedure are compared with 
the standard saddle-point approximation (see previous section) in Fig. || for m = 0.001 (small mass) and n, = n = 1 
({n)=N). 

We may now ask if a levelling in the free energy occurs for small masses m. As we wrote in section 4.1 this takes 
place for n < nxtop- However, due to the screening of the topological charge by the fermion determinant for small 
masses we expect \top ~ 777 (see Appendix C) so that the levelling is ruled out for reasonable values of n. We note 
that the macroscopic limit is reached only for Nra large. For N = 200 (Nm < 1) many modes are still missing in 
the sum ([[5]) (dotted line). For N = 10 4 (Nm 3> 1) and 9 < n/2, the quasi-analytical procedure (solid) and the 
saddle-point approximation (dashed) are in agreement. 

Finally, we now ask whether the standard saddle point method fails for 9 > tt/2, leading possibly to a cusp at 
9 < 7r, contrary to expectations. For Nf = 1 we can calculate the fixed asymmetry partition function exactly (|l^) and 
use it for a comparison with the quasi-analytical procedure (see above). The outcome confirms the standard saddle 



'Since (n) is fixed even, the periodicity is tv and not 2tv as noted in section 3.1. 
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point result with no cusp at 9 = it, and infirms the quasi-analytical result for large values of 9. We conclude that the 
saddle-point approximation carried prior to the x-resummation is only valid for 9 small (the 1/N terms at large 9 are 
important), while the one carried after the x-resummation is valid whathever 9. 

We have numerically checked, that most of the present observations carry to Nf > 1. In particular, a cusp may form 
in the latter for sufficiently degenerate quark masses, in agreement with the saddle point analysis discussed above. 



V. CONCLUSIONS 



We have analyzed the effects of a finite vacuum angle 9 on the vacuum partition function described by a matrix 
model, both in the quenched and unquenched approximation. The results are subtle. 

In the quenched case, we have found that the free energy exhibits a cusp at finite 9 that is sensitive to the precision 
of the numerical analysis. On this point, we are in agreement with the lattice analysis g. However, we have further 
noticed that the results are also sensitive to the maximum density of winding modes n. For a small compressibility or 
a peaked distribution in n, n is similar to n*, the mean winding density. For large enough n > 1 the position of the 
cusp is moved to 9 = n for high enough precision. This observation may be of relevance to the lattice results |t],|^] . In 
this context, it would be interesting to compare the lattice distributions for n± in @,| to the gaussian ones we have 
used in our work. 

In the unquenched case, a similar dependence on n is found, where n is also interpreted as the maximum density 
of zero modes. For sufficiently large n and large masses the quenched results are recovered. Each mass decouples at 
to ~ x„, although for large and degenerate masses to ~ NfX* ( m units where the quark condensate is one). For Nf = 1 
and small masses with mV > 1, the screening of the topological charge takes place, and the saddle point solution 
holds without any cusp. For Nf — 2, 3, the numerical results are found to agree with a saddle-point analysis, and 
results from anomalous Ward identities and effective Lagrangians. A cusp at 9 = n occurs for sufficiently degenerate 
quark masses. We have found that the use of the saddle-point approximation at large 9 requires care. 
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APPENDIX A: ALTERNATIVE SADDLE-POINT WITH N F = 1 

An alternative saddle-point analysis can be directly performed for fixed \ using representation (^). For that, we 
define 

X = -in-y. (Al) 

The solution for the saddle point equations gives 



—to 2 ± -^/to 4 + 4(to 2 — y 2 ) iy 



2z z 



, -to 2 ± ^to 4 + 4(to 2 -y 2 ) iy 

t-sp - Yz 7 • {A2) 

For a peaked distribution in n we expect n ~ (n). The P sp 's are related to the condensate by 

1 „ . „ n 



M = -d m log Z = y (Pj p e ie + P sp e- W ) . (A3) 



The value of y is fixed by requiring the vanishing of the term proportional to x- Hence, the consistency condition 
reads 

log P+M + 2i0 + = or arctan |- + 9 + ^ = (A4) 

D-iy x* D x* 



7 



where D = (to 2 ± -^/to 4 + 4(to 2 — y 2 ))/2, with y satisfying 

(n+ -n-}=^y (i* e* - P 5p e- <e ) = n* y. (A5) 

In the consistency condition, the principal branch of the logarithm is retained, making the saddle-point result mani- 
festly 2tt periodic in 0. We note that the present derivation is equivalent to performing the saddle-point calculation 
for y without substituting the form flA2| ) for Pg P 's. The associated free energy is 



F=-UogZ(0) = ~^[hr A 



D 2 +y 2 



l + y 2 2 n* 2 

-m y . (A6) 



D 2 +y 2 x* 



The subtracted free energy is VAF(9) = -In Z(0)/Z(O). 



APPENDIX B: TOPOLOGICAL DENSITY 

The topological density (n+ — n_) measures the difference between the number of zero modes with plus and minus 
charges, in the volume V fixed by the width of the quenched topological susceptibility x*- This is also the amount of 
U(l) charge in the vacuum state thanks to (|T^). At the saddle point (Nf = 1) 

(n+-n_) = id e \nZ{9) /V = {Pi p e l9 ~ P sp e- W ) =n»». (Bl) 

We see that y measures directly the topological density at the saddle point. For a large mass to 

2/=-X-^n -XJU-^] (B2) 

while for a small mass to 

"v — 2n 

± y = -to tan 01 cos 9 1 - to 2 — sin 20 + 0(m 3 ) . (B3) 

The ± solutions correspond to the transformation 0^—0, whereas the absolute value corresponds to the change of 
the branch in the solution of the consistency equation. 



APPENDIX C: TOPOLOGICAL SUSCEPTIBILITY 

The topological susceptibility measures the variance of (n + — n_) in the vacuum state at finite 0. It is simply 
Xtop = —dy/d9 n*. In the large mass limit 

Xtop = X* ( 1 - J + °( m ^) > ( C1 ) 

giving xtop = X* in the quenched case (m = oo). In the small mass limit and |0| < 7r/2, 

Xtop = n»TO I cos 01 + ii*to 2 — cos 20 + C(to 3 ) . (C2) 

2x* 

For > 7r/2 we have — n--) 2 ) < 0, which is possible since the measure is not semi-definite. 
The quark condensate for large masses is 

_ i{q 1 q) = 2i U _ ^0 2 ) (C3) 
to \ iijto/ / 

and for small masses is 

( 1 x* — 2n* \ 

tm/t\ n oc\^Q — n m I — -I- — - — Qin 2 f) 1 -I- (Km 2 } 



i(qlq) = n* cos0 - n*TO - + — sin 2 + C(to 2 ) (C4) 

2 2x* 
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The anomalous U(l) Ward identity is given by 



Xtop = -im(q f q) - m 2 (q^ 5 q q^ 5 q) 
For small masses, the insertion of (|C^C4|) into (|c|) yield 



(C5) 



Xps 



Xps — n * \ 2 



211, 



2X* 



■ cos 2 9 



(C6) 



for the pseudoscalar correlator (rightmost term in (C5)). It is finite in the chiral limit. For 9 = we recover the 
result |25f] . For large enough n* > x* a t certain finite angles the pseudoscalar correlator becomes zero. For large 
masses we have \ps = (n* — X*)l m2 - The present relations generalize readily to Nf > I. 



APPENDIX D: RESOLVENT FOR FIXED X 

The resolvent for the matrix model considered here reads 



G "^(^ (iJ w - > 



ime~ ie 



Since the overlap matrix elements do not mix different flavors the resolvent splits into a sum of 1-flavor resolvents 
and we get effectively Nf copies of the appropriate 1-flavor eigenvalue distributions. 

Since the matrix is nonhcrmitian it turns out that the eigenvalues lie on a curve (more precisely on two intervals - 
see below). This comes from the decomposition (Nf = 1) 



.<///< W \ . „ , , / — msin^ W 

m sin 



Wt ime-* 8 I = tmcos0 1+ I wl ™«^ti ) (° 2 ) 



So the eigenvalues are just the eigenvalues of the two-level hermitian chiral system displaced by imcos8. We will 
write z' — z — im cos 6 and introduce the self energies defined by 



1 \ ( 7-V 



ime lB W \ I \ 
W' ime -1 



l 



(D3) 



Taking into account the fact that the random matrices are asymmetric, we obtain the following equation for the 
self-energies 

5 1 = ~ 1 ~\ v (D4) 

z' — msmu — 1j2 

1 + x 

5 2 = —± (D5) 

with x — x/N- The trace of the resolvent is just (Si + E2V2. Hence 



n( , , 1 V 1 t + U msm8 

G(z) = z ^— x (D6) 

where a — z' 2 — m 2 sin 2 9. The eigenvalues lie on two intervals determined by 

z' 2 ut -m 2 sin 2 0= X " (D7) 

2 ± V4 - x 2 

and, for asymmetric matrices, there are additional Dirac delta spikes at z = m sin 9 when x is negative and at 
z = — msin^ when x is positive. Explicitly the eigenvalue distribution reads (A = Xji + 
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u(X) = S(Xi - tocos6>)<! x s {Xb, + msin0) + ^-(Xr - msin0W 1- - + ^\ (D8) 

Z7T V O (7 



for x positive and 



u(X) = S(Xi -m cos 6) {\x\5(X R - msinfl) + — (X R - msm9)\ 1 - - + %\ (D9) 
^ Z7r V a cr z ) 

for x negative. The structure of these distributions is not universal, but may be useful for understanding the 9 
structure from the bulk QCD spectrum using cooled lattice gauge configurations. 
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FIG. 2. Free energy from a CP" model (left) § and (right) §. 
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